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Alkaline  fuel  cells  are  low  temperature  fuel  cells  for  which  stationary  applications,  e.g.  cogeneration  in 
buildings,  are  a  promising  market.  In  order  to  guarantee  a  long  life,  water  and  thermal  management 
has  to  be  done  in  a  careful  way.  In  order  to  better  understand  the  water,  alkali  and  thermal  flows,  a  two- 
dimensional  model  for  an  Alkaline  Fuel  Cell  is  developed  using  a  control  volume  approach.  In  each  volume 
the  electrochemical  reactions  together  with  the  mass  and  energy  balance  are  solved.  The  model  is  created 
in  Aspen  Custom  Modeller,  the  development  environment  of  Aspen  Plus,  where  special  attention  is  given 
to  the  physical  flow  of  hydrogen,  water  and  air  in  the  system.  In  this  way  the  developed  component,  the 
AFC-cell,  can  be  built  into  stack  configurations  to  understand  its  effect  on  the  overall  performance.  The 
model  is  validated  by  experimental  data  from  measured  performance  by  VITO  with  their  Cell  Voltage 
Monitor  at  a  test  case,  where  the  AFC-unit  is  used  as  a  cogeneration  unit. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Fuel  cell  technology  permits  to  generate  emission-free  electric¬ 
ity,  using  hydrogen  as  fuel  and  producing  water  and  heat  instead  of 
polluting  exhaust  gases.  The  alkaline  fuel  cell  (AFC)  was  the  first  to 
be  put  in  practical  service  in  space  applications.  However,  despite 
of  its  early  success  the  current  research  and  development  activities 
are  more  directed  towards  proton  exchange  membrane  fuel  cells 
(PEMFC)  as  the  technology  of  choice  for  vehicular  applications  [1  ]. 
Also  more  research  is  devoted  to  high  temperature  fuel  cells  as  the 
solid  oxide  fuel  cell  (SOFC)  and  molten  carbonate  fuel  cell  (MCFC) 
which  show  some  interesting  prospectives  in  new  power  cycles 
for  electricity  plants  [2-4].  In  a  changing  climate  of  energy  supply 
and  demand,  smart  grids  with  distributed  generation  as,  e.g.  micro¬ 
cogeneration  can  be  the  future  for  our  energy  market  [5].  Compared 
to  other  micro-cogeneration  technology  fuel  cells  show  the  high¬ 
est  electrical  efficiency,  and  are  promising  systems  for  residential 
and  small-scale  cogeneration  applications  [6].  The  main  reason  for 
its  high  performance  is  its  modular  construction:  therefore  its  effi¬ 
ciency  is  hardly  influenced  by  its  size.  [7]  In  spite  of  the  limited 
development  efforts  for  AFC  technology  it  is  still  competitive  with 
other  fuel  cells  for  power  applications.  Further  improvement  in  the 
AFC  technology  is  still  possible  and  will  strengthen  this  position 
[1].  Intensys* 1 *  has  developed  an  AFC-system  which  can  be  used 


*  Corresponding  author  at:  KHKempen,  University  College,  Kleinhoefstraat  4, 
2440  Geel,  Belgium.  Tel.:  +32  473  53  92  25;  fax:  +32  14  58  48  59. 

E-mail  address:  ivan.verhaert@khk.be  (I.  Verhaert). 

1  Intensys  BVBA,  Hoge  Mauw,  Arendonk,  Belgium,  http://www.intensys.com. 
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as  a  cogeneration  unit.  With  help  from  the  VITO  the  system  is  CE 
approved  and  commercially  available.  [8]  The  system  is  operational 
in  Vilvoorde  and  Hasselt  as  a  demonstration  project.  Experimen¬ 
tal  data  is  generated  to  examine  its  performance  in  time  and  at 
different  operation  modes.  Advancements,  regarding  overall  ther¬ 
modynamic  behaviour  and  overall  efficiency,  are  still  possible.  In 
this  work  a  model  of  an  alkaline  fuel  cell  is  built  in  Aspen  Custom 
Modeller.  This  model  focuses  on  the  thermal  management  of  the 
AFC  with  circulating  electrolyte.  With  this  model  several  system 
designs  can  be  examined  in  Aspen  Plus,  in  order  to  initiate  further 
improvement. 

2.  Thermodynamics  of  an  alkaline  fuel  cell 

2.1.  General  operation 

An  AFC  operates  by  introducing  hydrogen  at  the  anode  and  oxy¬ 
gen/air  at  the  cathode. 

•  At  the  hydrogen  inlet  a  gas  mixture  of  water  vapour  and  hydrogen 
enters  the  gas  chamber  of  the  fuel  cell.  The  hydrogen  diffuses  out 
of  the  gas  chamber  into  the  working  area  of  the  anode. 

•  At  the  oxygen  inlet  C02-free  air  or  pure  oxygen  arrives  in  the  gas 
chamber.  The  oxygen  diffuses  into  the  working  area  of  the  cathode 
to  take  part  in  the  reaction. 

Both  electrodes,  anode  and  cathode,  are  separated  by  a  circu¬ 
lating  electrolyte,  a  6M  potassium  hydroxide  solution  (Fig.  1).  At 
the  anode  hydrogen  reacts  with  hydroxyl  ions  into  water  and  free 
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Nomenclature 

hA  overall  conductance  (W  K-1 ) 

A  effective  area  of  the  fuel  cell  (m2) 

ENernst  thermodynamic  potential  (V) 

F  Faraday’s  constant  (C  mol-1 ) 

G  Gibbs  free  energy  (J  kg-1 ) 

h  total  enthalpy  (J  kg-1 ) 

I  load  current  (A) 

m  mass  flow  rate  of  species  i  (kg  s-1 ) 

Mi  molecular  weight  of  species  i  (kg  mol-1 ) 
pi  (partial)  pressure  (bar) 

P  power  (W) 

Q  heat  (W) 

R  universal  gas  constant  (J  mol-1  K-1 ) 

T  temperature  (°C) 

U  voltage  (V) 

v  velocity 

Re  Reynolds  number 

Nu  Nusselt  number 

q  coefficient  for  the  electrochemical  model 

Greek  symbols 

a  transfers  coefficient 

r\  overvoltage  (V) 

Subscripts 
an  anode 

cat  cathode 

FCB  fuel  cell  body 

e  electrical 

el  electrolyte 

surr  surroundings 

w  water/water  vapour 


electrons  (Eq.  (1)). 

H2  +  2(OH)~  -*  2H20  +  2e"  (1) 

Within  the  electrolyte,  the  water  is  transported  from  the  anode  to 
the  cathode.  An  external  electric  circuit  leads  the  electrons  to  the 


ALKALINE  FUEL  CELL 

Electrical  Current 


Fig.  1.  Working  Principle  of  an  Alkaline  Fuel  Cell  (source:  U.S.  Dept  of  Energy- EERE). 


cathode.  At  the  cathode  oxygen  reacts  with  water  and  electrons  in 
hydroxyl  ions  (Eq.  (2)) 

H20  +  |02  +  2e_  -*■  2(OH)“  (2) 

These  ions  flow  from  cathode  to  anode  through  the  electrolyte, 
to  sustain  the  total  electrochemical  reaction.  Combining  both  reac¬ 
tions  the  overall  reaction  (Eq.  (3))  shows  that  the  end  product  is 
water,  which  can  be  removed  in  one  or  both  gas  streams  or  in  the 
electrolyte,  depending  on  the  fuel  cell  configuration. 

H2  +  i02+2e-^H20  (3) 

The  overall  reaction  is  exo-energetic.  This  energy  has  an  electric 
part,  which  is  consumed  in  the  external  electric  circuit,  and  a  ther¬ 
mal  part,  which  results  in  a  temperature  rise  inside  the  fuel  cell. 
To  maintain  the  overall  fuel  cell  temperature,  heat  is  removed  by 
outlet  mass  flows  or  by  losses  to  the  environment. 

2.2.  Overview  of  AFC  models 

To  predict  fuel  cell  performance  in  order  to  improve  it,  mathe¬ 
matical  models  were  developed  to  predict  electrical  power.  In  the 
early  nineties  Kimble  and  White  proposed  a  model  for  a  complete 
fuel  cell,  where  they  take  into  account  the  polarization  and  physical 
phenomena  going  on  in  the  solid,  liquid  and  gaseous  phase  of  both 
anode,  separator  and  cathode  regions,  assuming  a  macro  homo¬ 
geneous,  three-phase  porous  electrode  structure  [9].  The  model 
divides  the  fuel  cell  in  five  layers,  a  gas  diffusion  layer  and  reaction 
layer  for  both  anode  and  cathode  and  a  separator,  containing  the 
electrolyte.  In  1999  this  model  was  the  basis  for  the  model  of  Jo  and 
Yi,  where  they  corrected  some  invalid  correlations  and  parameters, 
e.g.  in  the  open-circuit  potential  and  the  liquid  diffusion  parame¬ 
ters  [10].  Both  Kimble  and  White  and  Jo  and  Yi  used  immobilized 
or  re-circulating  electrolyte  and  removed  water  in  the  gas  streams. 
In  2006  Duerr  et  al.  translated  the  model  of  Jo  and  Yi  to  a  stack- 
model  in  a  Matlab/Simulink  environment  and  added  dynamics  to 
the  electrical  part  of  the  model  [11  ].  Few  details  were  given  on  the 
calculation  method  and  the  estimation  method  of  some  physical 
parameters.  These  models  are  all  meant  to  predict  the  polarization 
curve  or  electric  response  of  the  fuel  cell.  In  1998  Rowshanzamir 
et  al.  studied  the  mass  balance  and  water  management  in  the  AFC 
[12]  by  only  applying  mass  balances  and  diffusion  laws  (Stefan- 
Maxwell)  for  the  gas  diffusion  layer.  Flowever  their  model  does  not 
provide  any  prediction  on  fuel  cell  performance.  Because  of  the  fuel 
cell  type,  each  model  could  make  an  assumption  on  water  disposal. 
All  stipulated  that  water  was  removed  in  both  gas  streams  or  (for 
the  Orbiter)  in  the  hydrogen  stream.  The  fuel  cell  (system)  which  is 
subject  of  our  research  removed  water  by  both  exhaust  gases  and 
electrolyte  flow.  After  the  fuel  cell,  all  these  streams  are  collected  in 
one  reservoir,  where  eventually  at  nominal  working  point,  the  air 
flow  is  responsible  for  the  disposal  of  water  (vapour). 

2.3.  Model  development 

Using  the  same  methodology  of  H.  Fluisseune  et  al.  [13]  for 
PEMFC  a  lD/2D-model  with  control  volume  approach  is  obtained 
with  following  assumptions. 

•  Dynamic  pressure  losses  within  the  fuel  cell  are  neglected.  In  this 
way  the  total  pressure  can  be  assumed  constant  over  the  entire 
fuel  cell. 

•  The  temperature  is  assumed  to  be  uniform  in  each  control  volume 
and  all  output  flows  have  this  temperature. 

•  The  partial  pressures  within  the  control  volume  are  the  mean  of 
the  input  and  output  flow  in  the  second  dimension. 

•  The  heat  losses  from  the  gas  chambers  to  the  environment  are 
neglected,  because  the  heat  exchanging  surface  is  relatively  small. 
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All  heat  losses  to  the  environment  are  therefore  modelled  as  heat 
losses  of  the  fuel  cell  body  to  the  environment. 

The  three  parts  which  were  considered  in  the  control  volume 
model,  are  the  anode  and  cathode  gas  chambers  and  the  fuel  cell 
body,  where  the  reaction  takes  place  (Fig.  2).  The  model  is  modularly 
built.  In  this  way  a  more  detailed  model  can  be  obtained  by  serially 
connecting  several  individual  models.  Fig.  2  shows  a  screen  shot  of 
an  individual  model  in  Aspen  Custom  Modeller. 

For  each  control  volume  a  mass  and  energy  balance  is  used. 
The  energy  balances  are  closed  by  heat  fluxes  which  are  exchanged 
between  neighbouring  control  volumes  or  between  a  control  vol¬ 
ume  and  the  environment.  This  heat  flux  is  modelled  as  a  convective 
heat  flux.  Next  to  heat  and  mass  transfer  between  control  volumes, 
also  electric  energy,  generated  in  the  fuel  cell  body  is  transferred 
to  the  environment.  For  every  control  volume  the  mass  and  energy 
balances  are  specified.  Next  an  electrochemical  model  is  used  to 
describe  the  electric  behaviour  of  the  fuel  cell. 


vapour  from  the  cathode.  While  oxygen  diffuses  to  the  fuel  cell  body, 
water  vapour  diffuses  into  the  gas  stream. 

^O  2,  Fuel-in,  cat  —  ^O  2,  Oxygen  +  ^02,  Fuel-out,  cat 

™w,  Fuel-in,  cat  +  W-w, vapour  =  ^w,  Fuel -out,  cat  (6) 

2, Fuel-in, cat  =  ^N2,  Fuel -out,  cat 

Energy  balance: 

^ i™i,in,cat  '  hi, in, cat  +  Qcat  =  ^il33i0UtCat  '  ty, out, cat 

Qcat  =  QFCB,cat  (7) 

QFCB,cat  =  hAFC:catXTFCB  ~  T'cat ) 

The  cathode  outlet  temperature  is  one  of  the  main  parameters  on 
which  the  model  is  validated. 


2. 3 A.  Anode  gas  chamber 
Mass  balance: 

The  fuel,  a  mixture  of  hydrogen  and  water  vapour,  enters  the  gas 
chamber,  where  the  hydrogen  diffuses  to  the  fuel  cell  body  to  react. 
A  part  of  the  water,  formed  during  the  reaction,  diffuses  as  water 
vapour  back  into  the  hydrogen  gas  chamber.  A  mixture  of  unused 
fuel  and  water  vapour  leaves  the  gas  chamber  to  a  next  stage. 


2.3.3.  Fuel  cell  body 
Mass  balance: 

Within  the  fuel  cell  body,  the  driving  electrochemical  reaction 
takes  place.  The  mass  and  molar  balance  relates  the  hydrogen  and 
oxygen  consumption  to  the  generation  of  water  and  electric  current. 
The  current  is  linked  to  the  molar  flows  by  Faraday’s  law. 


^H2,Fue/-i'n,an  2, Fuel-out, an  —  ^H2,  Hydrogen 

™w, Fuel-out, an  ~  ™w, Fuel-in, an  —  W-w, vapour 

Energy  balance: 

^ i™i,in,an  '  h i,in,an  ~  ^ i^i,out,an  '  hi, out, an  =  Qan 

~OFCB,an  =  Qan 

QFCB,an  =  hAFQ  an  -(TfcB  ~  Tan  ) 


(4) 


(5) 


™  ,  Mw  •  lcen  _ 

mw,  el-in,  FCB  +  ^  p  ~ 

mw, el-out, FCB  +  mw,cat  +  mw,  an 
™  _  Mw  •  I ceu 

niH  2, Hydrogen  —  ^  p 

Mw  •  Iceii 


2,  Oxygen  — 


4  F 


(8) 


As  boundary  condition  in  the  simulation  the  fuel  entering  the  (first) 
gas  chamber  is  pure  hydrogen  and  there  is  no  hydrogen  leaving 
the  (last)  gas  chamber  as  unused  fuel,  because  it  is  an  end-of-pipe 
system. 

2.3.2.  Cathode  gas  chamber 
Mass  balance: 

The  entering  air  contains  oxygen,  nitrogen  and  water  vapour. 
There  is  a  large  excess  of  air  because  air  is  used  to  remove  water 


Energy  balance: 

Within  the  fuel  cell  body  the  different  layers,  gas  diffusion,  cat¬ 
alytic  and  separator  layer  are  enclosed  [9-11].  Although  different 
layers  exist  in  the  fuel  cell,  in  this  paper  the  properties  of  the  elec¬ 
trolyte/separator  are  taken  to  define  the  thermodynamic  behaviour 
of  the  fuel  cell  body.  The  mass  flows  between  fuel  cell  body  and 
gas  chamber  however  will  only  consist  of  gas  in  accordance  to  the 
boundary  conditions  of  the  gas  diffusion  layers.  This  will  affect  the 
enthalpy  of  these  stream  and  will  limit  the  mass  flow,  because  the 
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partial  pressure  cannot  exceed  the  saturation  pressure. 

E imi,in,FCB  '  ^i,in,FCB  = 

£ imi,out,FCB  '  hi,out,FCB  +  QfCB  +  Pe 
QfCB  =  QFCB,sur  +  OfCB,cat  +  QFCB,an 


(9) 


QFCB,sur  =  hAFCB,sur-(TFCB  ~  Tsur *) 

The  energy  balance  of  the  fuel  cell  body  consists  not  only  of  incom¬ 
ing  and  outgoing  mass  streams  (with  certain  enthalpy)  and  heat 
flows,  but  also  of  an  electric  power  output.  This  output  is  more 
detailed  in  the  electrochemical  model. 


Table  1 

List  of  the  used  semi-empiric  parameters. 


2.3.4.  Electrochemical  model 

The  goal  of  the  presented  research  is  to  find  a  model  which  can 
predict  outlet  temperatures,  next  to  electrical  output.  Therefore  the 
electrochemical  model  is  based  on  a  parameter  approach  of  the 
polarization  curve,  which  relates  current  to  voltage  and  operating 
parameters  and  is  able  to  show  effects  of  composition,  flow  rate, 
temperature,. . .  on  the  cell  performance  [14].  The  electric  power  is 
a  product  of  both  current  and  voltage.  Current  can  be  determined 
using  Faraday’s  law.  The  fuel  cell  voltage  is  defined  by  the  Nernst 
potential,  the  activation  overvoltage,  the  ohmic  voltage  and  the 
diffusion  losses. 


Pe  =  UI  (10) 

Id  —  ^Nernst  ~  hact  ~  Pres  ~  Pdiff  (1*0 

The  Nernst  potential  is  calculated  out  of  the  Gibbs  free  energy  from 
the  electrochemical  reaction.  Gibbs  free  energy  is  calculated  using 
the  properties  and  libraries  of  Aspen  environment. 


r  _  A  Go  ,  RTcell 

Nernst  2f  +  2f 


Wph2)+  2ln(Po2) 


(12) 


The  activation,  ohmic  and  diffusion  losses  are  calculated  by  the 
equations  found  in  literature  [15,16].  To  calculate  these  losses  fol¬ 
lowing  expressions  ((13),  (14)  and  (15))  are  used. 


P  act  — 


RT 
a  n-F 


In 


(13) 


Tires  —  Re  '  I  cell 


(14) 


(,5) 

The  activation  overvoltage,  given  by  Eq.  (13),  is  a  function  of  the 
exchange  current,  j0,  which  is  temperature  dependent  given  by  Eq. 
(16). 

Jo  =  Ci  ■  exp  (  Jr— )  (16) 

V  1  cell  J 

The  ohmic  resistance  is  a  result  of  resistance  due  to  resistance  of 
electron  collector  and  a  term  inversely  related  to  conductivity  of 
the  electrolyte.  This  conductivity  is  related  to  the  cell  temperature. 
This  relation  results  in  a  linear  expression  (Eq.  (17)) 


Re  =  c3  -c4  ■  Tceu 


(17) 


2.3.5.  Discussion  on  parameter  estimation 

The  empiric  model  shows  great  similarity  to  the  one  used  by 
Amphlett  et  al.  [17]  and  later  by  Huisseune  et  al.  [13]  for  their  mod¬ 
elling  on  PEM  fuel  cells.  The  general  equations  are  useful  for  all  types 
of  fuel  cells.  The  parameters  however  will  differ.  All  physical  prop¬ 
erties  (enthalpy,  Gibbs  free  energy,. . .)  are  linked  to  the  libraries  of 
the  Aspen  environment.  A  number  of  parameters  however  are  semi- 
empiricly  obtained.  Table  1  shows  all  the  semi-empiric  parameters 
for  this  model.  In  this  paragraph  a  brief  discussion  on  every  param¬ 
eter  is  included. 

At  a  certain  current  density,  jL,  the  diffusion  losses  will  be  too 
high.  The  value  is  based  on  stack  construction  (See  Section  3.1), 
maximum  current  and  active  cell  area  (A  =  0.064m2).  Ref.  [16] 


°  0,0  0,2  0,4  0,6  0,8  1.0 

Air  Flow  Rate  (kmo/hr) 


Fig.  3.  Experimental  set-up. 
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shows  for  other  fuel  cell  types  values  of  20,  000  Am-2.  The  rela¬ 
tively  low  value  is  due  to  the  use  of  air  instead  of  oxygen  (reduction 
by  a  factor  5)  and  due  to  the  low  development  rate  in  AFC  tech¬ 
nology.  The  activation  parameters  are  a  best  fit  for  the  measured 
data,  within  an  acceptable  range.  Those  ranges  went  for  a  from  0.1 
to  0.5  and  for  j0  from  10-4Acm-2  or  lAm-2  to  10-9Acnrr2  or 
10-5  Am-2.  The  resistive  parameters  are  in  the  same  order  as  those 
in  [17,13]  for  PEMFC.  For  the  convective  values  we  can  refer  again 
to  Amphlett  et  al.  [17]  and  Fluisseune  et  al.  [13].  Similar  results  are 
shown  for  a  stack  of  the  same  size.  Both  anode  and  cathode  trans¬ 
fer  coefficient  are  in  a  similar  range.  The  transfer  coefficient  to  the 
surroundings  is  higher  because  here  the  system  consists  out  of  4 
stacks  and  in  [17,13]  only  one  stack  is  considered  which  is  more 
compact.  In  a  first  model  approach  the  convective  heat  transfer 
coefficient  was  fixed,  although  it  is  function  of  amongst  others, 
temperature  and  velocity.  Especially  this  last  influence  parameter 
seemed  inneglectible.  In  the  model  validation  it  is  shown  that  at 
high  air  flow  rates  (overflow)  the  model  underestimates  the  air  tem¬ 
perature.  Therefore  both  convection  coefficients  are  configured  as  a 
function  of  the  gas  flow  rate  (Eq.  (18)).  This  resulted  in  much  better 
experimental  validation  (See  Section  4). 

hAcat  =  hAan  =  C5  •  if6  (18) 

3.  Experimental  work 

3.1  Experimental  set-up 

Experimental  tests  for  model  validation  were  done  at  the  VITO 
on  the  AFC-system  of  the  KFILim.  Fig.  3  shows  a  schematic  view  of 
the  experimental  set-up  [8].  The  system  consists  of  4  stacks.  The 
stacks  are  supplied  with  ambient  air  at  a  measured  temperature 
(measuring  point  a).  The  air  flow  is  regulated  by  an  electric  fan 
with  variable  speed.  The  air  flow  rate  is  indirectly  measured  by  the 
electric  current  absorbed  by  the  fan  motor  (point  b).  Earlier  mea¬ 
surements  show  a  linear  correlation  between  this  current  and  the 
air  flow  rate  in  the  nominal  working  range.  The  air  is  guided  through 
a  C02 -scrubber,  where  a  slight  temperature  rise  is  expected  (point 
c  is  estimated).  The  air  is  parallel  divided  over  the  stacks.  Within 
every  stack  the  air  is  again  parallel  divided  over  the  different  cells. 
The  air  overflow  will  be  guided  to  the  electrolyte  tank,  but  because 
other  configurations  will  be  examined  this  is  not  yet  incorporated 
in  the  drawing.  The  outlet  air  temperature  is  measured  (point  d). 
A  centrifugal  pump  pumps  the  electrolyte  (point  e)  out  of  an  elec¬ 
trolyte  reservoir  and  enables  the  circulation  of  the  electrolyte.  The 
electrolyte  flow  is  indirectly  measured  (point  e)  in  the  same  way 
as  the  air  supply,  by  a  measured  correlation  between  flow  rate  and 
current  absorbed  by  the  electro  motor  of  the  pump.  Before  the  elec¬ 
trolyte  enters  the  stacks,  its  temperature  is  measured  (point  f).  Then 
the  electrolyte  is  in  parallel  divided  over  the  different  stacks  and 
cells.  After  the  passage  through  the  fuel  cells,  the  temperature  is 
measured  (point  g)  and  the  electrolyte  is  cooled  in  an  external  heat 
exchanger.  The  hydrogen  (99.98%  purity)  is  delivered  as  pressur¬ 
ized  gas  to  the  system  at  ambient  temperature  (point  a).  It  runs 
serially  through  the  stacks.  To  prevent  blockage  by  water  droplets, 
the  direction  will  be  alternated  (not  shown  in  drawing).  Finally  the 
electric  current  and  cell  voltages  are  measured  by  the  cell  voltage 
monitoring  (CVM)  (point  h)  [18].  First  every  cell  is  electrically  in 
parallel  connected  with  three  neighbouring  cells.  These  groups  of 
four  cells  are  than  serially  connected  in  a  stack  (24  groups  or  96 
cells  in  each  stack).  Finally  the  (four)  stacks  are  serially  connected 
to  each  other.  An  external  device  will  convert  the  DC-current  to  AC. 
To  validate  the  model  the  average  cell  and  total  system  values  were 
used.  All  the  described  measurements  are  logged  with  a  time  step 
of  one  second.  Several  steady  state  points  were  gathered  by  keep¬ 
ing  the  current  constant  and  averaging  KOFI-temperature  within 


Table  2 

List  of  all  maximum  standard  deviations  for  each  dataset. 


Parameter 

Gmax 

Total  Current 

2A 

Total  Voltage 

2  V 

Air  Flow  Rate 

0.183  kmolh-1 

Air  Outlet  Temperature 

2  °C 

Electrolyte  Flow  Rate 

2.2  kmol  h-1 

Electrolyte  Temperature  Inlet 

4°C 

Electrolyte  Temperature  Outlet 

4°C 

acceptable  range.  The  temperature  is  controlled  by  a  cooler  in  the 
secondary  circuit,  which  cools  the  electrolyte.  Not  every  working 
point  could  be  reached  because  the  capacity  of  the  cooler  was  lim¬ 
ited.  The  combination  of  high  currents  and  low  temperature  was 
therefore  not  possible.  Secondly  at  low  currents  it  was  not  pos¬ 
sible  to  reach  high  temperatures  as  there  was  no  external  heater 
in  the  electrolyte  circuit.  Because  the  goal  of  this  study  is  to  get 
a  model  for  normal  working  range,  these  limitations  were  accept¬ 
able. 

3.2.  Data  analysis 

All  gathered  data  -  which  were  obtained  during  1  week  of  exper¬ 
iments  -  were  chronologically  obtained  and  small  fluctuations  were 
inevitable.  To  reduce  the  amount  of  data  to  a  workable  set  of  work¬ 
ing  points,  all  data  were  stored  and  divided  in  several  datasets.  The 
boundary  condition  to  every  dataset  is  that  current  and  temper¬ 
ature  are  at  least  10  s  within  a  certain  range  (±2.5°  C  and  2.5  A) 
without  interruption.  For  each  dataset,  an  average  value  and  stan¬ 
dard  deviation  is  calculated  for  each  parameter.  If  the  standard 
deviations  for  the  measured  parameters  are  within  the  limitations 
as  described  in  Table  2,  the  dataset  is  considered  as  a  valuable  work¬ 
ing  point  for  model  validation.  This  led  to  a  50  valuable  working 
points,  containing  information  about  electrolyte  flow  rate,  inlet  and 
outlet  temperatures,  about  air  flow  rate  and  temperature,  about  cell 
current  and  voltage. 

3.3.  Discussion  on  experimental  results 

Voltage  and  current:  fuel  cell  performance 

To  examine  the  performance  of  the  fuel  cell  system,  all  experi¬ 
mental  data  are  shown  in  a  current-voltage  diagram  or  polarization 
curve  (Fig.  4).  The  horizontal  axis  shows  the  total  current  of  the  fuel 
cell  system,  which  is  a  measure  for  the  current  density.  The  vertical 
axis  shows  the  total  voltage  of  the  fuel  cell  system.  As  expected  the 
voltage  drops  with  rising  current. 

Fuel  cell  performance  and  the  electrolyte  temperature 

To  examine  the  influence  of  the  electrolyte  temperature  the  data 
were  divided  in  three  data  sets.  Each  data  set  represents  a  number 
of  data  points  within  a  certain  temperature  range  of  the  electrolyte 
input  flow. 

•  The  blue  dots  represent  experimental  data  with  an  electrolyte 
temperature  between  25 0  C  and  40  °C. 

•  The  red  squares  represent  data  at  a  temperature  between  40 0  C 
and  60  °C. 

•  The  green  triangles  represent  the  experiments  at  relatively  high 
temperature  from  60 0  C  to  75  °C. 

After  comparing  these  three  sets  of  data,  can  be  concluded  the 
voltage  is  increased  with  temperature  (Fig.  4). 
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Fig.  4.  Polarization  curve  for  different  electrolyte  temperature  ranges  (experimen¬ 
tal). 
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Fig.  6.  Air  flow  in  function  of  electric  current  at  different  temperatures. 


Production  of  useful  heat 

The  electrolyte  flow  is  kept  almost  constant.  The  rise  in  temper¬ 
ature  is  therefore  a  good  indication  of  the  useful  heat  production. 
In  Fig.  5  the  influence  of  temperature  and  electricity  load  at  useful 
heat  production  is  illustrated.  The  vertical  axis  shows  the  temper¬ 
ature  rise  of  the  electrolyte,  which  is  representative  for  the  heat 
production.  The  horizontal  axis  shows  the  total  electricity  produc¬ 
tion  (DC)  of  the  fuel  cell  system.  The  data  points  shown  on  the  graph 
are  divided  in  three  groups,  according  to  their  temperature  range. 
This  is  similar  to  the  discussion  on  the  polarization  curve. 

At  high  temperature  a  linear  relation  between  electricity  gener¬ 
ation  and  heat  production  is  noticeable.  At  lower  temperature  this 
linear  relation  is  not  clear,  although  it  would  be  expected.  At  lower 
temperature  the  small  changes  in  environment  temperature  are  of 
a  bigger  influence  than  at  higher  temperature.  Comparing  high  and 
low  temperature  the  rise  in  temperature  will  be  higher  at  low  tem¬ 
perature.  This  is  expected  because  less  heat  is  lost  in  air  flow  and 
losses  to  the  environment. 

The  air  temperature  and  flow 

The  air  flow  rate  fluctuates  very  much  and  the  data  points  them¬ 
selves  also  show  a  relatively  large  standard  deviation.  The  air  flow 
is  a  regulating  parameter  in  the  water  and  temperature  control  and 
therefore  could  not  be  fixed;  a  certain  fluctuation  was  inevitable. 


Electrical  Power,  Pc  (Watt) 


Comparing  air  flow  rate  and  total  electric  current  at  different  elec¬ 
trolyte  temperatures,  a  linear  correlation  is  shown  at  temperatures 
below  60  °C,  despite  the  large  fluctuations  (Fig.  6). 

At  higher  temperature  this  correlation  also  exists  but  the  flow 
rate  is  increased  for  better  cooling  and  removal  of  the  water  vapour. 
Above  and  under  60  °C,  the  designed  working  temperature,  a  dif¬ 
ferent  control  strategy  is  used. 

As  shown  in  Fig.  7  the  air  temperature  follows  the  electrolyte 
temperature.  No  influence  of  the  different  air  flow  rates  for  the  high 
temperatures  was  noticed.  (See  also  next  section) 

4.  Model  validation 

The  model  was  verified  with  experiments  on  stack  level.  Current, 
hydrogen  input  temperature,  electrolyte  input  flow  and  tempera¬ 
ture  and  air  input  flow  and  temperature  were  input  data  for  our 
model  and  experimental  set-up  (See  Section  4)  and  the  model  is  val¬ 
idated  by  comparing  the  fuel  cell  voltage  and  output  temperatures 
for  both  electrolyte  and  air.  A  number  of  modelled  and  measured 
data  points  are  shown  to  illustrate  the  results  in  Table  3.  The  dif¬ 
ference  between  the  original  and  improved  model  has  only  got  an 
impact  on  the  air  temperature.  Only  at  the  discussion  of  this  param¬ 
eter  the  models  will  be  compared  to  each  other.  First  the  model  is 
validated  at  voltage  and  electrolyte  temperature. 
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Fig.  5.  Useful  heat  (rise  in  electrolyte  temperature)  compared  to  power  output. 


Fig.  7.  Relation  between  air  and  electrolyte  temperature  (experimental). 
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Table  3 

Number  of  modelled  and  measured  data  points. 


Data  Points 

Experimental 

Model 

Improved  Model 

nr 

I stack 

A 

Tin, el 

°C 

U stack 

v 

Tel,  out 

°c 

T air, out 

°c 

U stack 

v 

T elec,  out 

°c 

T air,  out 

°C 

Tf stack 

Telec,out 

T air, out 

1 

18.5 

29.4 

75 

31 

28 

75 

31 

30 

75 

31 

28 

2 

28.3 

33.4 

73 

36 

32 

73 

36 

35 

73 

36 

31 

3 

33.3 

37.9 

72 

41 

36 

72 

40 

40 

72 

40 

35 

4 

39.7 

54.0 

71 

57 

47 

71 

55 

51 

71 

55 

48 

5 

43.4 

34.0 

70 

37 

31 

70 

38 

36 

70 

38 

33 

6 

48.2 

53.2 

71 

56 
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53 

70 
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7 

53.4 

39.4 

69 

42 

35 

68 

44 

38 

68 

44 
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8 

58.5 

49.5 

68 
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9 
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Fig.  8.  Polarization  curve  for  different  electrolyte  temperature  ranges  (simulation). 

The  voltage  and  the  electrolyte  temperature 
The  model  shows  the  same  trends  for  temperature  dependency: 
with  a  rise  in  temperature  a  rise  in  voltage  is  achieved.  The  simula¬ 
tion  results,  which  are  presented  in  Fig.  8  show  a  good  resemblance 
compared  to  the  experiments. 

The  air  temperature 

The  linear  correlation  between  electrolyte  and  air  temperature 
however  is  not  confirmed  in  the  original  model  (Fig.  9). 

At  low  temperature  an  overestimation  was  made  and  at  high 
temperature  (above  60  °C)  the  model  underestimates  air  temper¬ 
ature.  A  detailed  look  at  the  difference  between  modelled  and 
measured  air  temperature  shows  a  linear  correlation  between 
model  difference  and  air  flow  rate  (Fig.  10). 

Possible  explication  is  the  positive  effect  of  the  flow  rate  on  con¬ 
vective  heat  transfer  of  the  air  channels  and  cathode,  which  was 
not  taken  into  account  in  the  original  model.  The  improved  model 
uses  a  correlation  (Eq.  (18))  between  convective  heat  transfer  and 
velocity.  This  is  based  on  the  definitions  of  the  number  of  Reynold 
(Eq.  (19))  and  Nusselt  (Eq.  (20))  and  the  form  of  their  correlation 
(Eq.  (21 )).  This  form  is  most  found  in  literature  (e.g.  See  References 
in  [19]). 

Re~v  (19) 

Nu  ~  h convection  (20) 

Nu  =  a.Reb  (21) 

The  improved  model  shows  a  better  result  for  the  prediction  of  the 
air  temperature  (Fig.  9). 


Electrolyte  Temperature(°C) 


Fig.  9.  Relation  between  air  and  electrolyte  temperature  (simulation). 
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Fig.  10.  Simulation  error  on  air  temperature  prediction  in  function  of  air  flow  rate. 

5.  Conclusion 

The  model  shows  good  results  in  predicting  the  two  most  impor¬ 
tant  output  parameters,  electrolyte  temperature  and  power  output, 
for  integrating  the  AFC  into  a  micro-cogeneration  unit.  Improve¬ 
ments  are  made  in  the  prediction  of  air  temperature,  which  resulted 
in  an  acceptable  model.  The  model  designed  in  Aspen  Custom 
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Modeller  can  now  be  used  to  integrate  in  several  energy  systems, 
where  an  energy  and  exergy  analysis  of  the  system  can  be  per¬ 
formed. 
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